Examples of nonuniform limiting distributions for the quantum 

walk on even cycles 

Malgorzata Bednarska, 1 Andrzej Grudka, 2 Pawel 
Kurzyriski, 2 Tomasz Luczak, 1 and Antoni Wojcik 2 

1 Faculty of Mathematics and Computer Science, 
Adam Mickiewicz University, Umultowska 87, 61-614 Poznan, Poland. 
2 Faculty of Physics, Adam Mickiewicz University, 
Umultowska 85, 61-614 Poznan, Poland. 

Abstract 

In the note we show how the choice of the initial states can influence the evolution of time- 
averaged probability distribution of the quantum walk on even cycles. 
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The analysis of discrete quantum random walks initiated by Aharonov et al. and 

,□11 

its possible applications for constructing efficient quantum algorithms (|2j-jy|) has recently 
attracted a lot of attention. Although many questions in this area remain open, it is well 



known that the behaviour of classical and quantum wa' 



seen by studying spreading, mixing and hitting times 



ks can be very different, as it can be 



El, Q and Q 



or limiting distributions 



3). One of the differences between quantum and classical walks we explore in this note is that 
one can start a quantum walk not from a single occupied node, but from the superposition 
of many nodes. The influence of the initial state on the behaviour of a quantum walk 
was studied by Tragenna et al. E3- I n El we mentioned the possibility of generating 
highly nonuniform limiting distributions in a quantum walk on even cycles starting from 
superposition states. In this note we show that the initial conditions can affect the time 
evolution of the total variation distance of time-averaged probability distribution in a decisive 
way. 

We shall study a quantum random walk on an even cycle with d nodes, using a model 
proposed by Aharonov et al. In this setting the nodes of the cycle are represented by 
vectors \v), v = 0, 1, . . . , d — 1, which form an orthonormal basis of the Hilbert space Hy. An 
auxiliary two-dimensional Hilbert space Ha (coin space) is spanned by vectors \s), s — 0, 1. 
The initial state of the walk is a normalized vector 

l*o) = ^2lsv\s,v) = ^2n/sv\s)\v) (I) 

s,v s,v 

from the tensor product H = Ha ® Hy. In a single step of the walk the state changes 
according to the equation 

\% +1 ) = U\%), (2) 

where the operation U = S(H (g) I) first applies the Hadamard gate operator H = 
75 s'(~ l) ss 'l s )( s 'l t° the vector from Ha, and then shifts the state by the operator 

S = J2\s)(s\®\v + 2s-l(modd))(v\. (3) 

s,v 

The operator U has been studied in p|, where we prove that 

U\4>jk) = c jk \(f)jk), (4) 
where the eigenvalues Cjk are given by 

c jk = -^((-l)Vl + cos 2 (27rj7d) -isin(27rj/d)), (5) 



for k = 0,1, and j = 0,1, . . . ,d — 1, while the corresponding eigenvectors are 

\<t>jk) = (%ifc|0) + a jk b jk \l)) ®^uj^\v), (6) 

where u>d = e 2m ^ d , 

a jk = l/yfd{l + \b jk \*), (7) 
b jk = ^((-l)Vl + cos 2 (27rj7rf) - cos(27Tj/d)) • (8) 

The probability distribution on the nodes of the cycle after the first t steps of the walk is 
given by 

p t (v) = j2\(sMm 2 - (9) 

As was observed by Aharonov at el. for a fixed v, the probability pt(v) is 'quasi- 
periodic' as a function of t and thus, typically, it does not converge to a limit. Thus, instead 
of Pt(v), the authors of 1] propose to consider time-averaged probability distribution 

1 - 

and its limiting distribution 

ir(v) = lim p t (v). (11) 

In order to present the global properties of the walk let us also define the total variation 
distance 

A * = ^E|ft(«)-5|. ( 12 ) 

v=0 

which measures how far is time - averaged probability distribution from uniform distribution. 
At tends to limit which will be denoted by Aoo = lim^oo A t . 

Figs, n 121 an d El show the evolution of the total variation distance A t for the case 
of three different initial states: l^jj ), \^^) and l^o^), which can be written as a 
superposition of some 2d, 2, and 4, eigenvectors \4>j k ), respectively. Thus, I^q 2 ^) = 
12jk9j^d V03 \ < j > jk) 1S the state with a single occupied node vq where gj k = aj k (l + ib* k )/\/2; 
l^o / = 75 (I03,o) + I09,o}) is a superposition of two degenerate eigenvectors; finally 

l*0 4) ) = 75 (I03,O> + |09, O ) - |015,O) " |021,o))- 

In the case a quantum walk starts with l^jj ), one observe decaying of the total variation 
distance to the nonzero value A^ (Fig. [TJ. An analytic form of A^ can be found in a 
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FIG. 1: Time evolution of the total variation distance At for the initial state l^n ) (A = ^4). 
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(2) 

FIG. 2: Time evolution of the total variation distance At for the initial state \^f ) [d = 24). 
Diamonds - numerical simulations, line - analytical value of Eq. (|23|) . 



9] 



similar way as in |°| (we remark that there is a minor error in the equation (22) in [9]). 



0: 



Thus, we get 



ir(v) 



d 



(13) 
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FIG. 3: Time evolution of the total variation distance A$ for the initial state I^q 4 ^) {d = 24) 
Diamonds - numerical simulations, line - analytical value of Eqs. (|12j) and (|26|), 

where 



vg x _ 1 
l-(-z)^ Z 0x0 d' 



(l + (_l)d/2) 



(i.e. , £ = 1 when d/2 is even, and £ = if d/2 is odd), and 

s = s(v) = min (\v — v \, d — \v — v \) 



(14) 



(15) 



(16) 



denotes the distance between nodes vq and v, and s' = d/2 — s. If d ^> 1 we can write A 2 ^ 
in a simple form. When £ = 

A£0 = l/d, (17) 



while in the case of £ = 1 



A (2d) = 2 _ £ / _ log 2 d-l/2 
00 d d 2 V log 2 ^ 



(18) 



where z = 3 — 2y / 2- For the particular case presented at Fig. d ()18j) gives the value 
Ad» = 0.054. Hence, if we start with a single occupied node, the total variation distance 
decreases steadily in time and its limiting value tends to zero as the graph size d grows. Let 
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us present now two examples of walk for which the dynamics of the total variation distance 
is dramatically different. Fig. |21 pictures the evolution of a walk where A t does not 
change in time. It starts with the initial state of the form the form 

-J= (\4>m,k) + \<Pd/2-m,k)) , ( 19 ) 



for m = 3 and k = 0. The state described by (jl9j) for m — 0, . . . , m max as well as the state 

(\<fid/2+m,k) + \<fid-m,k)) , (20) 

for m — 1, . . . , m max , (mma = |_(^ ~ 2)/4j) consists of two degenerated eigenvectors. Since 
the evolution of the superposition of any number of degenerated eigenvectors leads only 
to the global phase changes so the dynamics of the probability distribution is frozen and 
7i(v) = pt{v) = Po(v). For the states given by (|T5j) and (|2Uj) the limiting distribution takes 
form 

/ \ 1 (— l) w sina , . , . 

«(v) = 7] + 2 sin (a (2v + 1)) , (21) 

where a = 2-Km/d. Thus 

Af } = ^ 2 ^|sin (a( 2, + l))|. (22) 
c(V 1 + cos a 

When m divides d/2 the summation can be easily perform leading to 

Af° = % 7= ■ 1 2 (1 - cos (2a ( V + 1))) , (23) 
« V 1 + cos^ a 

where 

"-L^-ij. (24) 

For m = 3 and d = 24 ^ gives 0.204. 

The last example of the time evolution, depicted at Fig. El is, perhaps, most intriguing. 
The changes of the total variation distance in this case resembles the motion of the damped 
harmonic oscillator with shifted equilibrium. Let us emphasize also that the limiting value 
of the total variation distance is much higher than the initial one Aq 4 "* . The initial state 
|0q 4 ' ) ) is of the kind 

1 



.-, \ir'm,fe) + \(fid/2-m,k) ~ \4>d/2+m,k,k) ~ \ < t>d-m,k)') ■ (25) 
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The time-averaged probability distribution for the initial states of the form given by (|25|) 
can be described as 

-f \ Af \ i ui >n(2y? mfc (t + 1)) 

p t {v) = A{v) + B(v) — ^ , (26) 

where A(v) = ir(v), B(v) = (po(v) — vr(t>)) / sin(2<^ m fc), ip m k is the phase of the eigenvalue 
Cmk {cmk = e tiPmk ) and tt(v) is given by pijl . Fig. El presents Aj> calculated with the use of 
(|T2*|) and (|2*^|) as well as the results of numerical calculation. 

In conclusion, we demonstrated how the initial conditions affects the dynamics of the 
quantum walk on cycle. We gave examples for three different kinds of behavior of the total 
variation distance between given distribution and uniform distribution: decaying, constant 
and damped oscillating. 
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